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Abstract 

We introduce and analyze a system of two coupled partial differential 
equations with external noise. The equations are constructed to model 
transitions of monovalent metallic nanowires with non-axisymmetric in- 
termediate or end states, but also have more general applicability. They 
provide a rare example of a system for which an exact solution of nonuni- 
form stationary states can be found. We find a transition in activation 
behavior as the interval length on which the fields are defined is varied. 
We discuss several applications to physical problems. 

1 Introduction 

A locally stable system subject to random perturbations will eventually be 
driven from its basin of attraction by a sufficiently large fluctuation. The case of 
weak noise is particularly important, and often arises in stochastic modeling of 
physical and non-physical systems [TMlOj. Noisy systems in which spatial vari- 
ation of some intrinsic property cannot be ignored arise in numerous contexts 
and are therefore of particular interest [TT]. Examples include micromagnetic 
domain reversal [12, 13 , pattern nucleation 11 Hi . transitions in hydrogen- 
bonded ferroelectrics [T7], dislocation motion across Peierls barriers [TS], and 
structural transitions in metallic nanowires [191120] . 

Noise-induced escape can occur through either classical or quantum pro- 
cesses. In classical systems, where escape involves activation over a barrier, 
the source of the noise is often, but not always, thermal in origin. The field- 
theoretic techniques for computing escape rates in these systems were developed 
by Langer (211122] , who considered the homogeneous nucleation of one phase in- 
side another. In a quantum system at sufficiently low temperature, escape occurs 
by tunneling through a barrier; although the process is physically different, the 
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mathematical formalism is similar to that of the classical case. The basic theory 
here was worked out by Coleman and Callan [23 25 in the context of quantum 
tunneling out of a "false vacuum" . A review of the basic approach in both cases 
can be found in [9l[26] . 

All of these consider the case of an infinite system. The equivalent problem 
on a finite spatial domain was first studied by Faris and Jona-Lasinio [271128] . 
who developed a large deviation theory of the leading-order exponential term in 
activated barrier crossing for the special case of Dirichlet boundary conditions. 
General techniques for computing the activation prefactor were developed by 
Forman [3D] and McKane and Tarlie [3D]. In general, computation of the sub- 
dominant behavior of activated processes (the rate prefactor, the distribution of 
exit points, and related quantities) requires knowledge of the transition state, 
which describes the system configuration at the col, or barrier 'top'. When this 
state is spatially nonuniform, as typically occurs at all but the shortest system 
sizes, it can typically only be found as the solution of a nonlinear differential 
equation. This usually requires numerical techniques; nontrivial problems where 
analytical solutions can be found are rare. 

When spatial variation can be ignored, the dynamical evolution of a noisy 
system can be modelled using a stochastic ordinary differential equation; when 
it can't, a stochastic partial differential equation is required. This can and does 
lead to new phenomena; one of these is a phase transition in activation behavior 
as some system parameter is varied. In classical systems confined to a finite 
domain this transition can occur as system size (or some other parameter such 
as external magnetic field in a magnetic system) changes [31,32 ; in quantum 
systems, an analogous quantum tunnelingoclassical activation transition occurs 
as temperature increases [33-38 . This phase transition has a profound effect 
on activation behavior, and can have physically observable consequences [20] . 

Most of the cases studied to date require only a single classical field to 
describe the spatial variation of the system. However, situations can arise in 
which two (or more) fields are necessary to model the system. These have 
received little attention to date; one notable exception is the analysis of Tarlie et 
al. [39] which studied phase slippage in conventional superconducting rings. In 
this case, the spatial and gauge symmetries imposed by the physics allowed 
an exact solution of the transition state to be found. More generally, such 
symmetries are absent, and it is of interest to find more general systems in 
which exact solutions can similarly be found. 

In this paper we propose and analyze one such model that may have wide 
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applicability; it can be regarded as a generalization of the system studied in |39] . 
The model was motivated by the growing body of research on metallic wires of 
several nanometer diameters and lengths of the order of tens of nanometers 
[47] , which represent the ultimate size limit of conductors and are of interest from 
the point of view of both fundamental physics and technological applications. 
However, as we discuss in the conclusions, the model and its solutions may apply 
as well to other problems of interest. We therefore confine our discussion in this 
paper to the mathematical features of the model and its solution. Its specific 
application to nanowires — which will necessarily involve additional modelling 
— will appear elsewhere [JS] ■ 

2 Monovalent Metalic Nanowires 

Nevertheless, to set the stage for introduction of the model, we briefly discuss the 
problem of stability and lifetimes of nanowires composed of metals from either 
the alkali or the noble metal groups. Because the stresses induced by surface 
tension at the nanometer lengthscale exceeds the Young's modulus, such wires 
are subject to deformation under plastic flow |49j . A purely classical wire would 
therefore be subject to breakup due to the Rayleigh instability, and this in fact 
has been observed for copper nanowires annealed between 400 and 600° C [50] . 
However, a nanowire is sufficiently small for quantum effects to also play a 
role, and indeed a quantum linear stability analysis 1!) ol 53 showed that at 
discrete values of the radius, the Rayleigh instability is suppressed. These radii 
correspond to conductance "magic numbers" that agree with those observed in 
experiments [4^1IM1I55] . 

The linear stability analysis, however, ignores thermal noise that can induce 
rare but large radius fluctuations leading to breakup. A self-consistent approach 
to determining lifetimes [IS1I5B], which modelled thermal fluctuations through 
a stochastic Ginzburg-Landau classical field theory, obtained quantitative es- 
timates of alkali nanowire lifetimes, in good agreement with experimentally 
inferred values [46 [ I54 [ [55], The theory restricted itself to perfectly cylindrical 
wires, so that a single classical field could be used to represent radius fluctua- 
tions along the length of the wire. 

To test the assumption of axial symmetry along the cylinder axis, Ur- 
ban et al. |57j performed a stability analysis of metal nanowires subject to non- 
axisymmetric perturbations. They were able to show that, at certain mean radii 



3 



and aspect ratios, Jahn- Teller deformations breaking cylindrical symmetry can 
be energetically favorable, leading to stable nanowires with non-axisymmetric 
cross sections. They predicted that a typical mechanically controllable break 
junction experiment should observe roughly 3/4 cylindrical wires and 1/4 non- 
cylindrical. 

The mathematical problem can be understood as follows. Urban et al. [58) 
considered non-axisymmetric deformations in the wire cross-section of the form 
cos(to0), i.e., having m-fold symmetry. The radius function describing the 
surface of a wire of cross-sectional area irp(z) 2 at position z along the wire 
length is then 



R(p, <f>) = p(z) f 1 - A m (z) 2 /2 + A ™M cos[m(0 - cf> m )} , (1) 

\ym m J 

where the sums run over the positive integers and the deformation parameters 
\ m {z) that represent deviations from axial symmetry are considered small. 

Urban et al. studied deformations with to < 6, but focused mostly on m = 2, 
i.e., quadrupolar deformations. Deformations of order higher than m = 6 cost 
more surface energy and are therefore less stable, m = 1 corresponds to a 
simple translation. One can then consider only deviations from axisymmetry of 
the form shown in Fig. [TJ 




Figure 1: Solid line is the deformed cross section at to = 2 
Using a linear stability analysis (which again ignores large thermal fluctu- 
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ations) they found several sequences of stable wires, some with considerable 
deviations from axial symmetry. Fig. [5] shows the phase diagram of linear sta- 
bility in the configuration space of the two deformation parameters: the Sharvin 
conductance G s and the coefficient of the quadrupole deformation A2. The x- 
axis gives a measure of the average radius parameter p, which is related to the 
square root of the Sharvin conductance |59j . 
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Figure 2: Linearly stable quadrupolar Na wires (thick lines) at low tempera- 
ture, shown in the deformation configuration space (p, A2), where p is related to 
the Sharvin conductance by ^JG s /Gq ~ kpp/2. The thin grey lines show the 
thresholds for the openings of new conductance channels. Thick lines along the 
x-axis (A2 = 0) denote regions of stability for purely cylindrical wires. (From 
Urban et al. [58].) 

In order to construct a comprehensive theory of nanowire lifetimes, then, it is 
necessary to extend the theory developed in |19[l56j to non-axisymmetric wires. 
This requires, as noted earlier, a consideration of a classical stochastic Ginzburg- 
Landau field theory with two fields, with one representing variation of the radius 
along the longitudinal axis and the other the departure from axisymmetry. 

3 The Model 

From the preceding discussion, it is evident that a minimal description of fluctu- 
ations in a non-axisymmetric wire requires two fields: the first, which we denote 
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4>\{z), is related to p(z) and characterizes radius fluctuations about some fixed 
average po; the second, which we denote <f>2(z), is simply \2(z) and characterizes 
deviations from axisymmetry. A wire with a perfectly cylindrical cross-section 
everywhere would have 02 (z) = 0. 

Although developing a theory of fluctuations of non-axisymmetric wires pro- 
vides a physical motivation for studying classical Ginzburg-Landau theories with 
two fields, we are also interested in the general problem of such field theories. In 
this paper, we therefore construct and study a model which applies not only to 
some (but not all) transitions among non-axisymmetric wires, but also to other 
problems as well. This will be discussed further in Sect. [5] 

We therefore consider on [—L/2,L/2] two classical fields 4>i(z,t), 4>2{z,t) 
subject to the potential 

fffa.fc) = -f <t>\ + \& - f 4 + \4 + \m (2) 

where /ii and p,2 are arbitrary positive constants such that p\ 7^/12, breaking 
rotational symmetry between the two fields. A contour map of the potential is 
given in Fig. [3] 
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Figure 3: Equipotential contours in U(4>i,4>2) for the fields <p\ and (f>2 when 
fj,i = 3 and fj,2 = 2. 

If the two fields are subject to the potential ©, have bending coefficients 
of Ki and K2 respectively (in the nanowire case these can be related to surface 
tension), and are subject to additive spatiotemporal white noise, then their 
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time evolution is governed by the pair of coupled stochastic partial differential 
equations 



4>i = K X <t>'{ + mi 0i ~4>l~ 0102 + e 1/2 6(^ t) 

02 = + ^202 ~4>l- ^102 + ^^^(z, t) , 



(3) 



where < £,-(zi, t 1 )^ j (z 2 , t 2 ) >= 8{zi - z 2 )5{ti - t 2 )Sij, i,j — 1,2. We will make 
the simplifying assumption here that Ki = k 2 = 1. When fi± — fi 2 , it can 
be shown that the above equations are equivalent to those studied by Tarlie 
et al. |39| to model fluctuations in superconducting rings. The breaking of 
symmetry between the fields leads to an entirely different behavior from what 
was found for that case. 

The zero- noise dynamics satisfy 0; = —5H/8<f)i, with the energy functional 

H = |^ 2 2 (i(0i(z)) 2 + i(0^(z)) 2 + C/(0 1 ,0 2 ))dz. (4) 

The metastable and saddle, or transition, states are time-independent solutions 
of the zero-noise equations, [9 satisfying the Euler-Lagrange equations 

0'/ = -M101 + 0? + 0102 



At nonzero temperature, thermal fluctuations can drive the system from one 
metastable state to another. Such a transition proceeds via a pathway of states 
that first goes uphill in energy from the starting configuration, passes through 
(or close to) a saddle configuration, and then proceeds downhill towards the 
nearest metastable state. The activation rate is given in the T — > limit by the 
Kramers formula [9] 



Here AE is the activation barrier, the difference in energy between the saddle 
and the starting metastable configuration, and Tq is the rate prefactor. 

The quantities AE and To depend on the details of the potential ^ , on the 
interval length L on which the fields are defined, and on the choice of boundary 
conditions at the endpoints z — — L/2 and z = L/2. It was shown in |60) that 
Neumann boundary conditions are appropriate for the nanowire problem, and 
we will use them here. However, the theory is easily extended to other types of 




M202 + 02 + 0102 



(5) 



r - r ex P (-A£;/T) . 
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boundary conditions [ST], although for periodic boundary conditions care must 
be taken to extract the zero mode when performing prefactor computations |30] . 

From here on, we choose without loss in generality /ii > fi2- In this case 
there are two uniform metastable states <pi s = ±-J]I\, <j)2 S = and two uniform 
transition states (j)\ u = 0, 4>2u — ±\f~pa, as can be seen in Fig. El In the following 
section, we will see that the uniform transition states are true saddles, and 
therefore relevant for escape, only when L < L c = —j===. At L c , a transition 
occurs, and above it the transition states are nonuniform. 



4 The Transition State 

We have found analytical solutions to Eqs. ([5]) that describe nonuniform saddle 
configurations (hereafter referred to as "instantons" , in keeping with the usual 
practice). For general fix > /12 they are: 



fti.mi 2 ) = ±V™V ( 2 Mi - M2) - rn((j,i - jU 2 )sn(VMi - M2 z\m) (7) 



4%*m( z ) = ±VA*2 - rnSjii - ^ 2 )dn(VMi - 1^2 z\m) (8) 

where sn(.|m)and dn(.|m) are the Jacobi elliptic functions with parameter m, 
whose periods are AK(m) and 2K (m) respectively, with K(m) the complete 
elliptic integral of the first kind 62 . Imposing Neumann boundary conditions 
yields a relation between L and to: 

r (9) 

The instanton states for \i\ = 3, /12 = 2 and intermediate to are shown in Fig.SJ 



As to — > 0, L approaches its minimum length L c — —j===. In this limit, 
uf = 0, 020* — ^y/T^i & nd the instantons reduce to the spatially uniform 



v 2,0 

saddle states. As to — > 1~, L — > 00, and the instanton states become 



b[ n f = ±v£T tanh(VMi - M2 z) (10) 
t 

"2,1 



- ±*j2(j, 2 - Mi sech(VA»i - M2 *) (H) 



In the nanowire case, there is a nice geometric interpretation of this particular 
version of the escape process, which will be discussed in Sect. |6l 

At low temperatures, the leading order asymptotic dependence (AE in the 
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(a) 



(b) 



Figure 5: The instanton states 4>"m f° r infinite length as m — >• 1 for /Ltx = 3, 
= 2. 
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Kramers formula) of the escape rate can be computed as the difference in energy 
between the saddle and stable configurations, and is plotted in Fig. [HI 
As L — » oo, AE approaches ~ mHMi + 2/i 2 ). This is simply the energy 
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Figure 6: The activation barrier AE for Neumann boundary condition (solid 
line) at fix = 3, fj,2 = 2. The dashed line indicates the crossover from the uniform 
transition states to the instanton transition states at L c . The dots represent 
the numerical results of AE beyond L c . 

of the domain walls of Figs. [5a] [5b] 

5 Rate Prefactor 

For an overdamped system driven by white noise, the rate prefactor Tq can be 
derived in principle 9 , although this is often difficult in practice. The usual 
procedure is to consider small perturbations 771, i] 2 about the metastable state: 
4>i = 4>i. s + ?/i and <f>2 — 4>2,s + ?72- Then to leading order 77 = —Asfj, if = 
(771,772), where A s is the linearized zero-noise dynamical operator at 0i )S , 4>2,s- 
Similarly, A„ is the linearized zero-noise dynamical operator at the transition 
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state ^i u,^ u . Then 



(12) 

where X u .i is the single negative eigenvalue of A u , corresponding to the direc- 
tion along which the optimal escape trajectory approaches the transition state. 
Here Eq. (fl"2|) differs from the usual formula for Tq [2][9] by a factor of 2 in the 
denominator because we have two saddle points distributed symmetrically be- 
tween the metastable states (see Fig. [3J such that the transition can take place 
via either. In the next two sections, we consider two interval length regimes, 
each with different saddle configurations. 



Tn — — 1 



det A, 



det A„ 



5.1 L < L c 

In this regime, both the metastable and transition states are spatially uniform, 
allowing for a straightforward computation of Tq. Assume the system begins 
at the metastable state (— y/m, 0), passes through the transition state (0, y 7 /^), 
and finishes at (^/m, 0). Linearization around the metastable state gives 

- W = -A ( m ) = (~ d * + 2Ml ° )( Vl ) (13) 

dt\ m J s \ m J \ o ~dl + {m-m))\m) 

and around the transition state 

i(:M:)--C*~ ( 'r« v + J(:)- ™ 

The spectrum corresponding to A s is 

2m rei=0,l,2... 

n 2 =0,1,2... 



Tii = 0,1,2. 
n 2 =0,1,2. 




11 



When L < L c , all the eigenvalues of A„ are positive, except for A" 1=0 = 
— {Hi — H2), indicating that this is indeed a saddle configuration. The eigen- 
function corresponding to A^ 1=0 , which is spatially uniform, is the direction in 
configuration space along which the optimal escape path approaches (0, ^/jui)- 
The fact that the lowest positive eigenvalue =1 — > as L — >• L~ indicates a 
transition in the escape dynamics at L c . 

This in turn affects the rate prefactor, which for L < L c is 




rr 1=0 (2^ 



1 + 



L 2 



n 2 =0 



(2/X2 



\ 



UZ=o ((Mi 



M2) + -jx 



M2j 



A*i \/sinh(- v /^i — ^L) v sinh(^/2/iiL) 



^sinh(y2^2£) Vl sin(VAii - M2^)| 



(Mi ~ M2) ■ 



(15) 



The prefactor diverges at L c = 
vanishing of AJJ 1=1 



The divergence arises from the 



w — (/zi — fJ.2), indicating an appearance of a transverse 



soft mode in the fluctuations about the transition state. 



5.2 L > L c 

When the transition is nonuniform, computation of the determinant ratio is less 
straightforward. Here the linearized operator A u about the transition state is 



A u = 



-dl - (mi - 3(« 2 - (0^) 2 ) W&St 



2##$# -d 2 z - (M2 - 3(^) 2 - (0^) 2 ) J 

(16) 

and the linearized equations become a pair of second order coupled nonlinear 
differential equations. 

In order to calculate det A s /detA„ we make use of a generalization, due 
to Forman [29 , of the Gel'fand-Yaglom technique [53], suitable for differential 
operators in a 2 x 2-matrix. The Forman method is readily extendible to higher 
dimensions [64], and its central result is that 

dctA s _ dct[M + NY s (£/2)] 

detA„ " det[M + NY u (L/2)] ' ( ' 

The 4x4 matrices Y s (z), Y u (z) in (fTTf are "fundamental" matrices [55] ■ A 
fundamental matrix Y(z) has the property that, for any solution ff(z) of the 
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homogeneous equation Aff = 0, 



Y(z) 



{m(-L/2)\ 

mi-L/2) 

v'i(-L/2) 
\ V ' 2 (-L/2).J 



(18) 



The matrices Y s (z), Y u (z) are then defined to be fundamental matrices of 
the differential equations 



(19) 
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o\ 
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V2 











1 


m 


dz 
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V° 


Mi - M2 





V 


W 



d_ 

dz 



1)2 

W 






j.inst \2 



V 



inst ./.inst 

1,771^2,771 

inst \ 2 (^\ust \2\ 



i o\ 

1 




20M^ -(/i2-3(<o 2 -(<CT 2 ) V 



(20) 

which are just the equivalent first order versions of A s and A u (for L > L c ). 
The matrices Y s (z), Y u {z) will be discussed further below. 

The 4x4 matrices M and N encode the boundary conditions 



M 



ffi-L/2) 
J'{-L/2) t 



N 



'n{L/2) \ = 



(21) 



For Neumann boundary conditions, M and N have the forms 



M = 



/o 1 

1 



\0 0J 



(22) 



13 



and 



N 



/o o o o\ 



10 

\0 I J 



(23) 



Thus, the determinant ratio of two infinite-dimensional matrices is reduced to 
the ratio of the determinants of two finite-dimensional matrices. 



The most difficult part of the strategy outlined in the previous paragraph 
is computation of the fundamental matrix. We proceed as follows. The ma- 
trix Y(z) can be expressed as the product of H(z)H _1 (— L/2), where H(z) is 
constructed from the four linearly independent solutions of the first order equa- 
tion (|19p for the metastable state, or (J2DJ) for the transition state. Suppose 



h 2 (z) 
h[(z) 
\h' 2 {z)) 



(hz(z)\ 
hi(z 
h' 3 (z 

\K(z)J 



fh 5 (z)\ 

h e (z 



(h 7 {z)\ 

h 8 (z) 
h' 7 (z) 

\K(z)J 



are the four independent solutions, each satisfying (using the transition state 
for specificity) 



d_ 

dz 



(hx\ 

h 2 

w 











-(M2-3(^) 2 -( 



O 2 ) 



1 o\ 




1 


h 2 





h[ 


Oy 


\h'J 


(24) 





or equivalcntly 




-dl - ( Ml - 3(<^) 2 
20i n : 



« 2 ) 



-d 2 z - (p 2 
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and similarly for the rest. Then H(z) is: 



H(z) 



(hi(z) 
h 2 (z) 
K(z) 

\h' 2 (z) 



h 3 (z) 
hi(z) 
h' 3 (z) 



h 5 (z) 
h e (z) 
K(z) 
K{z) 



h 7 (z)\ 

h 8 {z) 

h' 7 (z) 

K(z).J 



(26) 



While it's elementary to obtain the independent solutions of A s ?f s = at 
the metastable state, there's no systematic way of dealing with nonuniform 
transition states. In our problem, we note (cf. Q170) that it is sufficient to 
compute the fundamental matrix Y(z) at the boundary z = L/2. One can 
therefore numerically integrate the coupled differential equations forward from 
z = —L/2 using (j2"0|) with four independent initial values, as follows. 

Suppose that (rji(z), 772(2), T)[(z), r)' 2 (z)) is any solution satisfying ([20]) . a, /3, 7, S 
are four arbitrary constants, and H(z) is the matrix of the four linearly inde- 
pendent solution vectors. Then it follows that 



fm(z)\ 



H(z) 




7 

w 



(27) 



Inverting, we get 



P 
7 

W 



fm(z)\ 

m(z) 

V'i(z) 
\v' 2 (z)J 



(28) 



Since a, /3, 7, 5 are arbitrary constants, we can write without loss of generality 



P 
7 

W 



= n-\-L/2) 





-L/2)\ 




-L/2) 




-L/2) 


[v' 2 (- 


-L/2)) 



(29) 
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Using this to replace the constants in Eq. (|27p. we arrive at 



mi*) 



■.U{z)U-\-L/2) 



fm(-L/2)\ 

%(-i/2) 
v'i(-L/2) 
\i 2 {-L/2)J 



(30) 



fl 








o\ 





1 














1 





V° 








V 



which is simply (|18j). It is now clear that only the boundary values of H(z) at 
z = ±L/2 are needed to compute Y(L/2). For example, we can choose 



(31) 



for H(— L/2), and numerically integrate forward using (|2"U1) to get H(L/2). Since 
each column of the identity matrix (|31|) is independent from the others, the 
columns in H(L/2) arc also linearly independent. In this way, Y(L/2) is readily 
obtained. 

Fig. [7] shows the prefactor To vs. L. The divergence of the prefactor as 
L — > L c from either side is preserved in the two-field case; its physical meaning 
has been discussed in [35]. The critical exponent characterizing the divergence 
is of interest and can be readily computed. When L < L Cl it is easy to see 
from that 

r ~ (L-L c )-*. (32) 

For L > L c , the critical exponent can be computed numerically and is also 1/2, 
as shown in Fig. [5J 



6 Discussion 

We have introduced a stochastic Ginzburg-Landau model with two coupled 
fields, and have found explicit solutions for the nonuniform saddle states that 
govern noise-induced transitions from one stable state to another. This model 
has a phase transition as system size changes, similar to what has been seen in 
the single-field case [3T1I31IBB1|B7] . 

The action functional dU) was designed to model nanowire stability and decay 
when non-axisymmetric cross sections appear during the time evolution of the 
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Figure 7: Prefactor Tq vs. L for /ij = 4, ^2 = 2. 




7: 



-79 -78 -77 -76 -75 

Figure 8: log To vs. log(i — L c ) near L+ for /zi = 4, /i2 = 2. The dots indicate 
numerical data. The line fit gives logTo = — 0.51og(L — L c ) + 1.35 for L very 
close to (and above) L c , 
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wire, either as beginning, ending, or intermediate states. However, the model — 
or its modifications — has wider applicability, which we briefly discuss below. 
We note also that the model can be viewed as a generalization of that used 
by Tarlie et al. [39] to model transitions between different conducting states in 
superconducting rings. 

As noted in Sect. [TJ the system analyzed here can serve to model certain 
transitions among different nanowire states, but not all. Nevertheless, it pro- 
vides a nice physical interpretation of the escape process in non-axisymmetric 
nanowires. The escape process occurs, as usual, via nucleation of a "droplet" 
of one metastable configuration in the background of the other. Fig. [5] indi- 
cates that the infinite nanowire begins as a cylindrically symmetric wire at one 
(uniform) radius and ends as a cylindrically symmetric wire at a different ra- 
dius, while passing through a sequence of nonaxisymmetric configurations. It is 
easy to see that if we were to take instead ^2 > Mii the model would describe 
a process where the starting and ending states have the same average radius 
but are non-axisymmetric with different deformation states, while the saddle is 
cylindrically symmetric. By varying the model parameters, one can interpolate 
between these two extreme cases. An extensive analysis of such transitions, 
including applications of the current model, will appear in |48j . 

The model presented here, or ones close to it, apply also to a variety of other 
situations. One of these is constrained dynamics of the sort that plays a role 
in viscous, or glassy, liquids (see, for example, j68][69]). The simplest situation 
corresponding to this sort of dynamics is described in Munoz et al. 70 . Fig. 2 
of that paper shows a particle attempting to diffuse from one stable position to 
another; if a second particle is in one of its two stable positions, it blocks the 
first. Aside from its natural occurrence in glassy liquids, such a situation may 
also be created and studied using an optical trap with metastable wells, such 
as that employed by McCann et al. [71] or Seol et al. [72] to test the Kramers 
transition rate formula under varying conditions. 

If one removes the diffusion terms from Eqs. ([3]), treating (f>i and (f>2 simply 
as spatial coordinates, then one has a model where two degrees of freedom act 
on each other in a similar manner, although in a mutual fashion (see, for exam- 
ple, |73|). However, inclusion of the diffusion terms describes a more interesting 
situation. Here one can think of the field 4>i(z) as describing a density of a 
certain particle species along an interval, and 02 {z) as representing a second 
species density. Diffusion of particles into or out of the interval depends on 
the state of 4>2 (z); its vanishing (with the set of parameters used in this paper) 
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"locks" the density of 4>i(z) to one of its stable configurations. In order for 
any diffusion to take place, §i{z) must change also, but its density distribution 
is coupled to that of 4>i{z). It would be of interest to pursue this further to 
investigate spatially inhomogeneous models of constrained dynamics. 

There are a number of other physical situations that can be stochastically 
modelled using two coupled classical fields. Examples of such models are the 
Fitzhugh-Nagumo model [711175] describing excitable media in general (and neu- 
rons in particular), and the Gray-Scott model [76] of chemical reaction-diffusion 
systems. These, however, are nonpotential systems, and can exhibit phenom- 
ena, such as limit cycles, that potential systems cannot. Nevertheless, if one 
examines transitions between "fixed points" of such systems, the leading order 
asymptotics should behave similarly to the model described here. On the other 
hand, the subdominant asymptotics, for example the prefactor, would require a 
different approach. 
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